Linear-optics realization of channels for single-photon multimode qudits 
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We propose and theoretically study a method for the stochastic realization of arbitrary quantum 
channels on multimode single-photon qudits. In order for our method to be undemanding in its 
implementation, we restrict our analysis to linear-optical techniques, vacuum ancillary states and 
non-adaptive schemes, but we allow for random switching between different optical networks. With 
our method it is possible to deterministically implement random-unitary channels and to stochasti- 
cally implement general, non-unital channels. We provide an expression for the optimal probability 
of success of our scheme and calculate this quantity for specific examples like the qubit amplitude- 
damping channel. The success probability is shown to be related to the entanglement properties of 
the Choi-Jamiolkowski state isomorphic to the channel. 



I. INTRODUCTION 

The most general transformation a quantum state can 
undergo is described by a quantum channel. For exam- 
ple, it may correspond to a controlled manipulation of 
a quantum system for some final aim — like in quantum 
information-processing protocols [l| — or it may represent 
an unwanted interaction with the environment. While in 
the first case implementing the respective quantum chan- 
nel is of direct practical interest, in the second case one 
may still be interested in the implementation of the chan- 
nel for the sake of understanding the role of noise, and 
how to counteract it, in real- world implementations of 
quantum information-processing protocols. It is worth 
remarking that striking effects in quantum information 
processing (QIP), e.g., the super-activation of the quan- 
tum capacity of channels 0, involve non-trivial noisy 
channels. 

Quantum optics is one of the best established phys- 
ical architectures for QIP 0-0] ■ I* has the advantage 
that the carriers of information — photons — interact nat- 
urally weakly with the environment, so that real noise 
is low. This makes simulating noise possible in a very 
controlled way. The workhorse of optical experiments 
is the manipulation via linear optical elements, such as 
beam splitters and phase-shifters. The linear-optics re- 
alization of channels has appeared in a number of works 
for specific cases. For example, random-unitary channels 
are common in experiments on decoherence-free and uni- 
tarily recoverable subspaces d, ^ and in the realization 
of mixed states [Tol . [Tlj . The simplest non-trivial exam- 
ple of a channel that is not random-unitary is perhaps 
given by the qubit amplitude-damping channel [l|. The 
counting statistics of this channel have been simulated 
using linear optics [T2l|, and a stochastic linear-optical 
implementation with a fixed success probability of 50%, 
independent of the value of the damping parameter, has 
been suggested [isj . 

In this article, we propose a linear-optics scheme for the 
stochastic exact realization of an arbitrary channel for 
single-photon multimode qudits. Under constraints moti- 



vated by the ease of experimental realization, our scheme 
achieves an optimal probability of success. An interesting 
result is that such a success probability is related to the 
entanglement properties of the Choi-Jamiolkowski state 
isomorphic to the channel [3, ■ This connection al- 
lows us to apply results in entanglement theory (l6| to 
the quite different problem of channel realization. 

Our results provide an optimal strategy for the realiza- 
tion of arbitrary channels, an important building block 
in experimental studies of QIP. In the specific case of 
the qubit amplitude-damping channel, our scheme pro- 
vides a significantly higher efficiency than alternative 
schemes JJJ without leaving the subspace of the encod- 
ing of the input state. In contrast to [I4I, this allows us 
to further process the output of the channel. 

The paper is structured as follows. In Section [TTl we 
provide definitions, fixing both the framework and the 
notation. In Section [llll we illustrate in detail the prob- 
lem we consider, that is, the realization of a quantum 
channel with a fixed set of tools. In Section IIV[ we pro- 
vide a scheme to realize any channel perfectly albeit only 
stochastically. In Section|Vl we relate the optimal success 
probability of the method proposed to the entanglement 
properties of the Choi-Jamiolkowski state isomorphic to 
the channel of interest. In Section IVIl we use this re- 
lation to provide bounds on the probability of success, 
both in the specific case of qubits, for which we are able 
to give analytic bounds, and qudits. In Section IVIIl we 
apply our technique to two examples, one being the qubit 
amplitude damping channel. Finally, we conclude and 
discuss possible future venues to investigate. 



II. DEFINITIONS AND FRAMEWORK 

The state of a quantum system may change over time 
due to some internal dynamics, to an interaction with its 
environment or to a measurement performed on it by an 
observer. Any physical transformation a quantum sys- 
tem can experience can be modeled as a quantum chan- 
nel A : Pin I— Pout- Every channel acting on a system 
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S admits a dilation, that means, it can be realized as 
some unitary interaction with an ancilla E, which is sub- 
sequently discarded: 

A[ps] = TlEiUsE PS «) CJE uIe), 

with oe the initial state of the ancilla More ab- 
stractly a quantum channel can be defined as a com- 
pletely positive trace-preserving linear map. Each chan- 
nel can be represented in the form A[p] = X^i^iP^I' 
where {Ai\ is a set of Kraus operators fulfilling the trace- 
preserving condition, ^l^i = I- The Kraus represen- 
tation of a channel is not unique. For instance, if {Ai} 
forms a Kraus decomposition of a channel A, the rela- 
tion Bi = UijAj, assuming Uij are the elements of a 
unitary matrix, will define a new decomposition {Bi} for 

We will frequently find the notion of operator norm 
useful in our discussions of quantum channels. Since we 
will always work in finite dimensions, the operator norm 
II A||oo of ^ corresponds to the largest singular value of A. 
An operator is an admissible Kraus operator — that is, it 
can be considered as part of some valid Kraus-operator 
set — as long as ||A||oo < 1- Any set of linear operators 
that satisfy the completion relation ^l^i = ^ ^^^^ 

constitute a valid quantum channel. 

In this paper, we will be interested in optical quantum 
systems. Each mode of an optical system is associated 
to a basis of Fock states \n), where n = 0, 1, 2... denotes 
the number of photons in the mode. The creation and 
annihilation operators, and a, respectively, provide 
a convenient notational framework for describing Fock 
states because of the relations a\n) — y/n\n— 1), a''\n) — 
\Jn -I- l|n-M), so that \ri) = (a^)" /V«T|0). These opera- 
tors have commutation relations [ai,aj] — Sij, [a|,aj] — 
and [oi, aj] = 0, where the indices i and j denote the op- 
tical mode and Sij is the Kronecker delta. 

The most common optical elements that are used in 
experiments for the manipulation of optical modes are 
beam splitters and phase shifters. Optical networks that 
are composed only of instances of these two elements are 
referred to as passive linear devices. Linear (quantum) 
optics is the part of quantum optics that, apart from the 
initial generation of entangled photon pairs and single- 
photon detection, deals only with passive linear devices 
[3|. Any unitary transformation U acting on d optical 
modes and preserving the total photon number can con- 
veniently be described by the way it transforms the cre- 
ation operators of the modes: 

outt \ ^ int /-|\ 

a, = 2^Uy a^ ', (1) 

where the elements of a unitary matrix. A trans- 

formation can be realized by linear optics if and only if 
it is of this kind. 

A phase shifter is an optical element that acts on a 
single mode as Ua^W^ = e^'^d^ . A beam splitter acts on 



two optical modes at a time and can be described by 

, ( cose -e'* sin 9 \ , . 

^""'^^ ^ \e-'^sine cos 9 )■ ^ > 

Any unitary that acts on d modes preserves the total 
photon number if and only if it can be implemented using 
these two devices [T7| . 



III. THE PROBLEM 

A qudit can be encoded by using one photon in d op- 
tical modes. An arbitrary logical state can be written as 
W = Z]f=i V'^Nl), with a logical basis {\iL}}f=i, where 
|*l) = o||0). We call this kind of encoding a d-rail encod- 
ing. This encoding is convenient when the interactions 
are limited to linear optics, because any unitary opera- 
tion ([T]) can be performed on the creation operators using 
linear optics, and under this encoding the basis states of 
a single qudit and the creation operators transform iden- 
tically. 

We are interested in the simulation of an arbitrary 
quantum channel A that acts on a qudit, using only pas- 
sive linear optics. What we want is a realization of A 
on the d-rail qudit, such that the logical subspace — the 
encoding — is mapped onto itself. This allows for further 
processing of the output of the channel. 

We will refer to the channel to be realized as the logi- 
cal channel^ to distinguish it from physical channels that 
evolve the state of the modes without necessarily pre- 
serving the logical subspace. 

As we noted earlier, we can always represent a channel 
in the form of a dilation where the channel is realized 
via the unitary interaction of the system with ancillary 
modes. We will limit ourselves to linear-optics evolution. 
For the sake of the ease of experimental implementation, 
we will assume several other reasonable restrictions: (i) 
to limit the number of photons that need to be generated, 
we only introduce ancillary modes that are initially in the 
vacuum state; (ii) in order to prevent the necessity of us- 
ing expensive feed-forward mechanisms (Pockels cells and 
high-speed high- voltage switches — see, e.g., 0,[I1]), we 
do not allow adaptive schemes; (iii) we will restrict our- 
selves to photon-number measurements, although it will 
actually turn out that commonly used threshold detec- 
tors suffice. 

When we consider the dilation representation of a 
channel, we can imagine that the final trace over the 
ancillary space corresponds to a measurement of the an- 
cilla, whose result is discarded. If we assume that the 
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ancilla starts in the (vacuum) state lO)^;, then we have 

Mps] 

= TTEiUsEPS^CTEUlE) 



d rails 



= E ( {j\e^Jm^Use\0) ) PS ( {j\e^Jm^Use\0) 
jk 



with {Mk}, Mk > 0, J2kMk = / a POVM on the an- 
cilla system E, and {\j)} an eigenbasis for aE- With 
our constraints — vacuum input ancillas and linear optics 
evolution — measuring the vacuum on the output ancillas 
is the only result that leaves the system within the en- 
coding. This can be seen easily considering the action 
of the linear optics unitary Use = C^lo on initial states 
|*l)|0)e, i = 1, . . . , d. We will consider d + e modes, with 
the first d used for the encoding, and the remaining e 
constituting the ancilla system E. Then we have: 



C/LoNL)s|0)E = C/Loal|0)s|0)£ 

d+e 

= ^u.,a]|0)s|0)B 



(d+e 
J=d+l 

From this expression it is evident that if we perform 
a photon-number measurement on the output ancillary 
modes and we obtain a result different from the vacuum, 
then the encoding is lost. The reason for this is that lin- 
ear optics preserves the photon number and the initial 
state of the system |iL)|0)E only has one photon in it. 
If the photon is measured in the ancilla, then the ini- 
tial state I^l) will be mapped out of the encoding to the 
vacuum, independently of which output ancilla mode the 
photon is measured in. 

Therefore, under the constraints that we have imposed, 
the only logical channels that can be realized determinis- 
tically must have a single Kraus operator. Such channels 
are necessarily unitary transformations, as it can be seen 
by the trace-preservation condition A = I. 



IV. THE SOLUTION: STOCHASTIC 
IMPLEMENTATION 

In this section, we will first see that any single logical 
Kraus operator (i.e., any Kraus operator of the logical 
channel) can be realized stochastically. Later we will in- 
troduce a further resource, randomness, and the ability to 
switch — according to such randomness — among different 
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FIG. 1: The diagram describes an optical circuit for a channel 
that would realize the Kraus operator A = VSU as well as a 
Kraus operator that would map the encoding to the vacuum. 
The boxes represent optical arrays that perform the unitary 
that labels them. The S transformation then consists of a 
set of beamspUtters, one for each mode, whose transmission 
coefficients are matched to the singular values of the matrix 
S. 



optical networks, and we will show that then any logical 
channel can be realized, albeit only stochastically. 



A. Implementation of a logical Kraus operator 

For any logical Kraus operator A that we want to apply 
to the input state, it is possible to construct an optical 
network such that A will correspond to the transforma- 
tion of the logical state if the output ancillary modes are 
detected to be in the vacuum state, given that they were 
in the vacuum state before the channel. Every Kraus 
operator has a singular value decomposition A = VSU, 
where U and V are unitaries and the matrix S is positive 
and diagonal, with diagonal elements < < I that 
correspond to the singular values of A. As unitary rota- 
tions can be realized deterministically on the encoding, 
in order to prove that A can be realized under our con- 
straints, it is sufficient to prove that any diagonal matrix 
S can be realized (see Figure [T]) . 

This is proven to be possible by considering the action 
([2]) of a beamsplitter on two modes. If the first mode, 
with creation operator a\ belongs to the encoding and 
the second mode is an ancilla — which means it starts in 
the vacuum — then the transformation that results when 
the vacuum is measured on the ancilla state effectively 
realizes the mapping i~j> cos 6'a'''. Since the angle 9 is 
arbitrary, we can simply implement any diagonal logi- 
cal Kraus operator S by using d ancillary modes and d 
beamsplitter, choosing the angles 9i such that Si = cos 6';. 



B. Perfect but stochastic implementation of an 
arbitrary logical channel 

A logical channel A that we may want to apply on 
the encoding will in general have a Kraus decomposition 
with n > 1. Therefore, by using a fixed linear 
optical network in the framework defined in Section Hill it 
will not be possible in general to simulate the channel, as 
only one logical Kraus operator can be realized per fixed 
optical network. 
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We will circumvent this problem by realizing individ- 
ually the various Kraus operators A^, i = l,...,n, in 
this way being able to preserve the encoding for each Ai . 
Roughly speaking, by randomly applying the different 
Kraus operators the logical channel A will be realized. Of 
course, this is possible only by allowing the linear optical 
network to change. We will introduce the possibility of 
switching among various optical networks — one for each 
Ai — according to a probability distribution {pi}- Each 
fixed optical network that we will introduce to realize 
the Kraus operator Ai will itself correspond to a quan- 
tum channel (see Figured]). This "average realization" 
of the logical channel will anyway be stochastic, because 
in the implementation of any Ai that is not unitary there 
will necessarily be a finite probability of ending up out- 
side the encoding, which corresponds to finding the input 
photon in the output ancillary modes. 

One important point is that, given the additional de- 
gree of freedom due to the choice of the probability dis- 
tribution {pi}, it is possible to consider the realization 
of a rescaled version Ai of Ai rather than exactly Ai. 
Of course each Ai must be a valid Kraus operator, i.e., 
Il^illoo < 1- We will use this rescaling degree of freedom 
to maximize the success probability for the realization of 
the channel. 

If we postselect on finding the output ancillary modes 
in the vacuum state, and if we choose the probabil- 
ity distribution {pi} and the Ai operators such that 
V^j^i = ^/PsuccAi for all i and for some < Psucc < 1, 
then the logical input state p will be mapped into the 
(unnormalized) logical state 

PiAipAl = psucc ^ AipAl . 



This will happen with probability Tr(^jPiJLip^J) — 
Psucc, and thus the logical channel A will be stochasti- 
cally implemented with probability Psucc (independent of 
the input p). 

Given that we want the channel to be realized perfectly, 
the figure of merit we care about is the probability of 
success Psucc, which we want to be maximal. One possible 
choice for the distribution {pi} and the operators Ai is 
trivially Pi = 1/n and Ai ~ Ai] this choice leads to 
a probability of success Psucc — ^/n. This strategy is 
independent of the properties of the Kraus operator {Ai} 
for the particular channel A, and depends only on the 
number of Kraus operators. As such, one can expect it 
to be non-optimal, and it certainly is in the case of a 
random-unitary channel 



K[p]=Y,q,U,pU}, 



with {Ui} unitaries and {qi} a probability distribution. 
Indeed, in this case an obvious better choice — and actu- 
ally optimal — is pi — qi, Ai = Ui, for all i, such that 




P„,cc A[p] in the 
encoding 

-fc^Outof the 

encoding 



FIG. 2: Pictorial representation of our scheme. Each solid 
rectangle represents a channel. The most external box is a 
mixture of the n channels F,; inside of it. Each of these in- 
ner channels corresponds to a linear-optics setup and for our 
scope its action on the encoding can be completely described 
without loss of generality by two Kraus operators, Ai and 
Bi. Each Ai preserves the d-rail encoding, while the Kraus 
operators Bi map an encoded state out of the encoding. If 
the condition .JplAi 



•^PsMccAi, for all i, is met, the overall 
result of randomly switching among the channels Fi accord- 
ing to the probability distribution {pi } is that of realizing the 
target logical channel A with probability Psucc (independent 
of the input p). 



The following theorem provides the optimal choice of 
the probability distribution {pi} and of the operators 
AiS to maximize Psucc, for any fixed Kraus decompo- 
sition {Ai}. 

Theorem 1. Given the Kraus decomposition {Ai} for 
the channel A, the optimal probability of .success for its 
realization is 



{{Ai}) 



This can be achieved by the choice pi = ^ ||!4°|j^ 
A., = II /ii Ai, for all i. 



(3) 



and 



Proof. From the condition ^/ptAi — ^/psucc Ai, for all i, 
one finds K > = Psucdl^JLi where we used the 

fact that 1 1 Ail loo < 1, because each Ai must be a proper 
Kraus operator. Summing over i and using X^iPi — 1; 
one arrives at Psucc < 1/ J2i 1 1^*1 loo- The probability dis- 
tribution and Kraus operators in the statement of the 
theorem saturate the inequality. □ 

Thus, the maximal probability of simulating the chan- 
nel adopting the Kraus decomposition {Ai} in our 
scheme is the inverse of J^i ll^illL- This quantity will 
in general depend on the specific Kraus decomposition. 
By optimizing over all Kraus decompositions we have the 
following. 

Corollary 1. (Optimal probability of success) In our 
scheme, the optimal probability of .success in the imple- 
mentation of A is 



Psucc(A) 



1 



{^.} Ei 



(4) 



Ps 



1. 



where the maximization is over all Kraus decompositions 
{Ai} of the channel A. 
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For convenience in the analysis to follow, we define the 
stochasticity of a channel as 



ct(A) = min > \\Ai 



(5) 



where the minimization is over all Kraus decompositions 
{Ai\ of the channel A, so that 

ftucc(A) = 

The name "stochasticity" is justified by the fact that the 
larger cr(A), the lower the probability of a successful re- 
alization of the channel. 

We remark that any specific Kraus decomposition will 
give an upper (lower) bound on the stochasticity (optimal 
probability of success) . 



V. RELATION WITH ENTANGLEMENT 
MEASURES 

The optimal success probability Psucc(A) of implement- 
ing a channel is clearly just a property of the channel it- 
self. Therefore it appears natural to look for a represen- 
tation of the channel that is independent of any specific 
Kraus decomposition. This can be done by considering 
the Choi-Jamiolkowski isomorphism. 

The latter is a one-to-one mapping between maps and 
operators [13, The isomorphism — explicitly in the 
direction from maps to operators — is defined as 



J(A) = (A®/) [^+] 

^\jZm^3\)®\i)^3\. 



(6) 



for some fixed choice of a maximally entangled state 

-J 

f 



(7) 



For our purpose, the interesting observation is that 
pure ensemble decompositions of the Choi- 

Jamiolkowski state Ja isomorphic to a channel A are in 
one-to-one correspondence with Kraus decompositions of 
A. This follows from the the fact that for any vector 
l^/j) e (g) there is an operator such that 

\^)^{A^®mi). (8) 

Here, the bar in \'\\)) denotes that the vector need not 
be normalized. In general, unless it is specified to the 
contrary with use of the bar notation, all states are 
assumed to be normalized. Thus, 



J(A) = ^p,|V.)(V'^| 



(9) 



E 



for \^,) = ^IV'^) = [A^^ ® /)|Vd )■ 

As we have seen, without the use of randomness the 
only channels that can be realized deterministically are 
unitaries, and using randomness and switching among 
optical networks we can extend this result only to random 
unitaries. Thus, we have that the only channels that 
can be realized deterministically in our framework are 
those whose Choi-Jamiolkowski state admits an ensemble 
consisting only of maximally entangled states. 

One then expects that channels whose probability of 
realization is high admit Kraus decompositions that are 
close to random-unitary. In turn this would mean that 
their Choi-Jamiolkowski states admit ensemble decom- 
positions that are highly entangled. We will show that 
this intuition is correct. 

The relation (|5]) implies 



if we consider the Schmidt decomposition 

SiV^K)K), with > 0, Y^i^i = (^1^), and A„ 
maxi{Ai}. Thus, we find for the stochasticity 



(10) 



|2 



cr(A) = min V II A, 

i 

= d min y^PiAmax(^/'j) 

{pi,4>i} 



d 1 



max 



dll- 



(11a) 
(lib) 

(11c) 

(Ud) 



where we used (|10p to move from the minimization over 
Kraus decompositions for A to the minimization over en- 
semble decompositions for J (A). The quantity 

Eaiip) = l-Xm.xW - l-max|(a,/3|V')P, 

where the maximum is taken with respect to factor- 
ized pure states |q;,/3) = |a)|/3), is the geometric mea- 
sure of entanglement for a bipartite pure state [l^ . 
More generally, for a multipartite pure state, the geo- 
metric measure of entanglement is defined as Eq(tP) = 
1 — max^^^p |((/)scp|V')P, with (/)sep a fully separable state. 
In the bipartite case, it coincides with the entanglement 
measure E2, which was defined in [2^ as one of a whole 
family of entanglement measures. The geometric mea- 
sure of entanglement has received a good deal of attention 
[21I [2^ because of its intuitive — even in the multipartite 
case — geometric interpretation as maximal overlap of the 
state of interest with a fully separable state, and because 
of its connections to other well-known entanglement mea- 
sures, like relative entropy of entanglement [1^ HJl. In 
the bipartite qudit case we are interested in here, one 
sees immediately that 



< Eg{iI^) < 1 



1 

d' 



(12) 
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The lower bound is achieved for a factorized pure state, 
while the upper bound corresponds to a maximally en- 
tangled state like the one in Eq. d?]). 

The geometric measure of entanglement is extended to 
the mixed-state case by the usual convex-roof construc- 
tion [H: 



(13) 



where we use ^ to stress that the resulting quantity is 
convex on the set of mixed states. 

The standard convex-roof is defined in terms of the en- 
semble containing, on average, the minimum amount of 
entanglement as quantified, in this case, by the geomet- 
ric measure of entanglement for pure states. Eq. (jlip 
involves instead the ensemble containing on average the 
maximum amount of entanglement. This corresponds to 
the concave-xooi construction 



'g(p) = max X^P^^^l-iAi), 



(14) 



where we use ^ to stress that in this way we are defining 
a concave function on the set of mixed states. 

For the sake of comparison with quantities better 
known in literature, let us mention that in the same way 
in which the entanglement of formation 26] Ep(p^^) = 
mm^p.,^AnyJ:^P^S{pt), with = Tr^ ( |^^^) (^^^ | ) 
and S{cr) = — Tr(cr log2 cr) the von Neumann entropy of 
a state a, is the paradigmatic example for a convex roof 
construction, the entanglement of assistance (27| 



(15) 



is the paradigmatic example for a concave-roof construc- 
tion. 

From pT|) it follows that the stochasticity is given by 



a(A) = d(l-i?g(J(A)), 



(16) 



and, as a result, the relation between the probability 
of success Psucc(A) for our scheme to realize a channel 
A and the entanglement properties of the related Choi- 
Jamiolkowski state J(A) can be expressed as 



VI. BOUNDS 

The evaluation of the stochasticity ([5]) for a given 
chaimel is in general a non-trivial computational prob- 
lem. The connection with entanglement that was de- 
veloped in Section |Vl more precisely Eq. p6)) . shows 
that calculating the stochasticity is equivalent to evalu- 
ating i?g(J(A)). In principle, this requires to check for 
all possible ensemble decompositions of J'(A), although 
one can use convexity arguments to restrict the search 
to ensembles of pure states for a Choi-Jamiolkowski 
state of rank r, similarly to the case of entanglement of 
formation (2^ . In this section we will be able to provide 
analytic upper and lower bounds that do not require any 
search. 

Entanglement of assistance and other concave-roof 
constructions have not been studied as well as convex- 
roof constructions. This is due to the fact that they are 
not entanglement measures [13, H^. Nonetheless, they 
are of interest because, e.g., they capture some proper- 
ties of multipartite entanglement. For example, the en- 
tanglement of assistance quantifies the average amount 
of entanglement that two parties — Alice and Bob — can 
share thanks to a measurement of a third party who holds 
the purification of the state. Thus, we will be able to 
make use of some results already derived in literature, 
in particular in [l^ and [l^ , to provide upper and lower 
bounds for the stochasticity a and the probability of suc- 
cess Psncc • 

We first start by illustrating the range over which Psucc 
can vary, illustrating the best and worst cases. We then 
identify a simple bound based uniquely on the mathe- 
matical properties of the operator norm. As we will see, 
such a bound will turn out to be pretty useful in investi- 
gating the examples of Section IVIIl We then proceed to 
consider bounds based on the entanglement properties of 
the Choi-Jamiolkowski state isomorphic to the channel 
of interest. 



A. Best and worst cases 

Given that Eg — and therefore Eq — satisfies ((H]), it 
follows from P7|) that 



succ (A) < 1. 



(18) 



Psucc(A) 



d(i-£;S(JA))- 



(17) 



We remark that, because Eq is a concave function on 
states, the probability of success Psucc is a convex function 
on channels, i.e., 

((1 - g)Ai + qk2) < (1 - g)Psucc(Al) -I- gPsucc(A2), 

for < 9 < 1- Of course, this could be concluded directly 
from dl]). 



As has been pointed out earlier, the upper bound in 
()18p can only be achieved by random-unitary channels, 
whose Choi-Jamiolkowski states can be written as con- 
vex combinations of maximally entangled states. The 
lower bound corresponds to Eq{J{K)) = 0, i.e., to the 
case where no ensemble for J(A) contains any entan- 
gled state. Such an occurrence was considered in the 
context of the study of the entanglement of assistance 
in [i^l , where it was proved that any state p^^ with van- 
ishing entanglement of assistance must be of the form 
p^B ^ \a){a\ (g) p^ or = p^ ® Given that we 
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are not considering general bipartite states, but states 
that are isomorphic to channels via the isomorphism ^ , 
for the first inequality in (fT8|) to be saturated it must be 
J(A) = |a)(a| (g) I/d. The latter condition implies that 
the output of the channel is a pure state independent of 
the input, i.e., A[p] = Tr(p)|Q;)(Q!|. 

It may seem strange that the channel that is almost 
the most trivial theoretically is the one that is the most 
difficult to implement under our constraints. One can 
provide the following intuitive explanation. The output 
state must be independent of the input, but at the same 
time still be in the encoding. Thus, the output state must 
include the photon of the input encoding, because the 
ancillary modes are initially in the vacuum state. This 
can be accomplished in the following way. In the scheme 
proposed in Fig. [U a random rotation is first applied to 
the input. Subsequently, d—1 of the encoding modes are 
measured while the remaining one is transmitted — that 
is, the transitivity of d — 1 of the d beam splitters is set 
to 0, while the remaining one is set to 1. Upon finding 
the vacuum in the measured modes, we know that the 
photon is in the only unmeasured mode, i.e., in some 
known logical basis state of the encoding. Then we can 
rotate such a state to the desired output state. Given 
the random rotation of the input, the probability that 
this procedure succeeds is exactly 1/d, independent of 
the input. 

B. Triangle-inequality bound 

By using the triangle inequality, it is straightforward 
to derive an upper limit on the success probability. 

Observation 1. (Triangle-inequality bound) For any 
quantum channel A, 

Proof. If {Ai} is any Kraus decomposition for the channel 
A then we have for the stochasticity: 

a(A) = min^P.II^ 

= min^||A,4||oo 

* i 

> min||VA,4|U 

i 

= l|A(/)||oo, 

where the inequality is due to the triangle inequality, and 
the dependence on the choice of the Kraus decomposition 
is lost because J^i AiAl = A(/), for any Kraus decompo- 
sition of A. □ 

This bound proves that it is necessary for a channel 
to be unital in order for us to implement it deterministi- 
cally using our scheme, because only for a unital channel 



||A(/)|loo = 1- This is consistent with the already argued 
fact that under our scheme only random-unitary chan- 
nels can be deterministically implemented. The bound 
is easily evaluated, being independent of any particular 
Kraus decomposition. 

We remark that any choice of a specific Kraus decom- 
position provides a lower bound on the probability of 
success. If such a lower bound matches the upper bound 
in (jI9p . then the given decomposition is proven to be 
optimal. 

C. Bounds based on entanglement properties of 
the Choi-Jamiolkowski state 

Now we will move to bounds that exploit the con- 
nection we observed between the success probability of 
our scheme and the entanglement properties of the Choi- 
Jamiolkowski state J(A). 

1. Qubit channels 

We will first focus on the qubit case. Not surprisingly, 
this is the case where we can employ most results from 
entanglement theory. In particular, we will mostly be 
concerned with the one of the most common entangle- 
ment measures, known as concurrence [lO, HH. The con- 
currence of a pure two-qubit state can be expressed as 

where 

with the complex conjugation taken in the computational 
basis and ay = . The definition extends to den- 

sity matrices via the standard convex-roof construction: 

C^(p)= min 

{pi,-4'i} 

It is straightforward to check that for a pure state the 
relation 

^G^i^) = \(}- v/l-CW') (20) 

holds. In [l^ it was argued that and are related 
E"^ = \ {\- VI - CU(p)2^ . We will instead be in- 
terested in the connection between E^ and the concave- 
roof of the concurrence, 

= max Vp^c(V^^). (21) 

I 

The examples of Section I VIII will prove that the rela- 
tion (|20l) does not hold for the concave-roof version of 
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the two quantities. Nonetheless, in order to obtain easily 
computable bounds for Psucc , we will exploit the remark- 
able fact that there is a closed expression for both (p) 
and C^{p). For the former it reads [3l| 

C7U(p) = max{0,Ai-A2-A3-A4}, 



where the A^'s are the eigenvalues of \fpP\fp in de- 
creasing order, with the state p = Oy® Oyp* ay ® Uy. For 
C"^(p) instead it holds 29]: 



(22) 



with F{<t,t) — Tr y\/\/o'Ty/aj the fidelity between two 

states a and r. We start by providing the following 
lemma that relates C'~'{p) and Eq{p). 

Lemma 1. Given any state p of two qubits, the following 
inequalities hold: 



C^ip) 



Proof. See Appendix [Xj 



(23) 



□ 



By using the lemma together with the relation ([TT]) for 
d — 2, and ([^. we immediately obtain the following 
result: 

Theorem 2. (Concurrence bounds) // J(A) is the Choi- 
Jamiolkowski state isomorphic to the qubit channel A, 
then 



F(J(A),J(A)) 



> Psucc(A) 



> 



(24) 



l + Jl-F{J{A),J{A)y 



symbol having probability p and the remaining d — 1 
symbols being equally likely. It coincides with the bi- 
nary entropy for d — 2. We remark that hd{p) is a 
concave function of p, and is monotonically decreasing 
for p > 1/d. This means that the inverse function 
h^^ : [0,log2d] 1] is weh defined. 

We are now ready to state the theorem that links en- 
tanglement of assistance and probability of success. 

Theorem 3. (Entanglement-of-assistance bounds) For 
a given qudit channel A, the following inequalities hold: 



2-Ea(JA) 



> Psucc(A) > 



1 



dh-^\E^{JK)) 



(26) 



where Ei^ is the entanglement of assistance, J(A) is the 
Choi-Jamiolkowski state isomorphic to the channel, and 
a ( A) is the stochasticity of the channel. 



Proof. See Appendix IB] 



□ 



VII. EXAMPLES 

In this section, we consider two examples for the qubit 
case: (i) the amplitude-damping channel and (ii) the 
probabilistic constant-output map. For these examples 
we are able to find analytic results for the probability 
of success, and we compare these exact results with the 
bounds we obtained in Section IVI CI The analytic re- 
sults are obtained by using the triangle-inequality bound 
of Observation [Hand the already remarked fact that any 
specific Kraus decomposition provides a upper (lower) 
bound on the stochasticity (optimal probability of suc- 
cess). 



2. Qudit channels and entanglement of assistance 

In the previous section we focused on the concurrence 
because its concave-roof version can be easily evalu- 
ated. Concurrence was generalized to higher-dimensional 
systems in a number of different ways |32h34| , and even 
high-dimensional "assisted" versions — i.e., concave-roof 
constructions — were considered [s^. As we mentioned, 
the most studied example of concave-roof construction 
is the entanglement of assistance (|15p . For this reason, 
we will provide bounds for the probability of success in 
terms of the entanglement of assistance. 

We will use the following generalization of the binary 
entropy that depends only on the number of possible out- 
comes, d, and one probability parameter, p: 



hd{p) ■■= -plog^p- (1 -p) log2 



1-p 
d-1 



(25) 



That is, hd{p) is the Shannon entropy of the probabil- 
ity distribution of d symbols (p, i ■ • • i )> with one 



A. Amplitude-damping channel 

The qubit amplitude-damping channel is used to model 
the decay of an excited state 1 1) into the ground state |0) . 
With probability e the channel causes the de-excitation 
of the input state. This de-excitation process is de- 
scribed by the Kraus operator Ai = y^|0)(1|. A sec- 
ond Kraus operator guarantees that the process preserves 
probability, i.e., that the channel is trace-preserving: 

^2 = |o)(o|-i-yr^|i)(i|. 

Using this specific decomposition and Observation [T] 
we find: 

||A(/)||oo<a(A)<||Ai||2^ + ||A2||L. 
For the lower bound, one finds 



|A(/)|U = \\AiA\ + A2Al\ 
1 



X 1 -e 
1 + e, 



e 
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FIG. 3: Optimal probability of success for the amplitude 
damping channel (color online): exact (solid line), and up- 
per and lower bounds in Eq. p7p (dashed and fine-dashed 
lines, respectively). 

and for the upper bound we get: 

ll^i||'oo + ll^2||L = 1 + 

Because these two bounds coincide, the Kraus decompo- 
sition is optimal, with a stochasticity of o'(A) = 1 + e 
and the optimal success probabihty Psucc(A) — (1 + e)^^. 
When e = the channel is trivially the identity channel 
and can be performed deterministically. However at the 
other extreme, e = 1, the channel becomes the constant 
map A(p) — Tr(p)|0)(0|, and it can only be realized with 
probability 50%. As we found in Section IVI A| the con- 
stant map is the map that has the lowest success proba- 
bility in our scheme. Therefore the parameter e that de- 
scribes the probability of de-excitation let us move from 
one extreme to the other of the stochasticity (or prob- 
ability of success). For the amplitude-damping channel 
we find that the concave-roof of concurrence (|2T|) satis- 
fies C"^(J(A)) = F(J(A), J(A)) = ^/T~e. We can then 
compare the analytic bounds from (|24p with the exact 
result we just found (see Figure [3]): 

\== > Psucc(A) = > (27) 

We remark that our scheme achieves a probability of 
success that depends on e and is close to 100% for e small. 
On the contrary, the scheme of [ll] has a 50% probability 
of success, independently of e. In our scheme, such a low 
success probability is just the worst case (e = 1). 

B. Probabilistic constant-output channel 

The second channel that we choose to analyze is a con- 
vex combination of the constant output channel and the 
identity map. Such a channel returns the input state 
with a probability of 1 — p or a fixed output state r with 
probability p. The map is then 

A : p i-> (1 — p)p + p Tt{p)t 

and its Choi-Jamiolkowski isomorphic state is simply 

J(A) = (l-p)|V2+)(V2+l+pr®^. 



FIG. 4: Optimal probability of success for the probabilis- 
tic constant-output channel, exact result and bounds (color 
online). Upper and lower surfaces (blue and green online, re- 
spectively) are the upper and lower concurrence bounds from 
Eq. (|24|l . evaluated for the Choi-Jamiolkowski state isomor- 
phic to the channel; the in-between surface (red online) is the 
exact success probability Psucc for realizing the channel. Be- 
cause Psucc is symmetric about the s — 1/2 plane we only plot 
this function for s G [1/2, 1]. 

We find the stochasticity of this channel by checking that 
the upper and lower bounds for the stochasticity that are 
generated from Observation 1 and a particular decompo- 
sition of the state match. Without loss of generality we 
can consider t to be diagonal in the computational basis, 
i.e., 

r = .|0)(0| + (l-.)|l)(l|. (28) 

In fact, only the degree of mixedness of r and not the 
specific basis influences the probability of success Psucc- 
This can be understood at the formal level by considering 
that a different choice of basis for r can be taken into 
account via a rotation U, which does not influence the 
entanglement properties of J(A) (see Appendix [C|) . 

Using this state, the bound stated in Observation [1] 
becomes: 

a(A) > ||A(/)|U = |l(l-p)/+2pT|U = l-p+2ps, (29) 

where we have assumed, without loss of generality, that 
s > 1/2. From this and by using Eq. we find 

£:g(J(A)) < 1/2 -p(s - 1/2). One can find an ensemble 
decomposition of J(A) that saturates the latter inequal- 
ity (see Appendix [P]) . therefore Psucc = (1 — P + 2ps)~^. 

This means that for this channel we also find that as 
the probability parameter p varies from to 1 we move 
from the identity map to a constant map. However, we 
can see from Figure 2] that the success probability of the 
constant map depends on how mixed the output state 
is. As expected from the discussion of Section fVI Al the 
lowest value for the success probability, Psucc(A) = 1/2, 
is only attained when the constant output state is pure. 

VIII. DISCUSSION 

We have provided a scheme to realize an arbitrary 
channel on a d-rail-encoded optical qudit, taking into ac- 
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count practical restrictions. In particular: (i) we only 
allow for operations that are realizable with high fidelity 
using linear optics; (ii) we only allow ancillary modes 
that are initially in the vacuum state, thus limiting the 
need for sources of single photons that are, as of now, 
still difficult to produce on demand; (iii) we do not allow 
feed- forward (i.e., adaptive schemes), which significantly 
reduces the cost of the necessary equipment and loss that 
is inevitably involved in such schemes due to the need of 
long fibers for optical delays; (iv) and we consider only 
photon-number measurements (actually, readily available 
threshold detectors suffice). The conventional linear op- 
tics toolbox (phase shifters and beamsplitters), as well 
as the possibility for randomly switching between differ- 
ent optical networks, are the only elements needed for 
the realization of our method. These restrictions ren- 
der our technique of immediate interest to linear-optical 
implementations that can be realized using state-of-the- 
art experimental techniques. Within this framework, it 
turns out that any channel can, in principle, be realized 
perfectly, albeit only stochastically. The only channels 
that can be realized deterministically are random-unitary 
channels. Given that post-selection is a commonly used 
technique in linear-optics experiments, this restriction ef- 
fectively only slightly reduces the success probability of 
an experimental realization, and we are able to provide 
an expression for the optimal probability of success. This 
probability turns out to be related to the entanglement 
properties of the Choi-Jamiolkowski state isomorphic to 
the channel of interest. More precisely, we were led to 
evaluate the "assisted version" of the geometric measure 
of entanglement, i.e., the concave-roof extension of the 
measure to mixed states, for such a state. While we are 
not aware of a closed formula for it, not even for two-qubit 
states, we were able to provide upper and lower bounds 
in terms of the concave-roof of concurrence (for qubits) 
and of entanglement of assistance (for general qudits). 

Besides tackling the problem of evaluating, in general, 
the concave-roof of the geometric measure of entangle- 
ment, i.e., the probability of successful realization of our 
scheme with the restrictions considered in this paper, fu- 
ture research will focus on the relaxation of said restric- 
tions, that is, on the analysis of more general schemes for 
the realization of channels. 

For example, the use of ancillary states that are not 



initially in the vacuum certainly improves the realization 
of certain channels. Indeed, we saw that the worst-case 
scenario is that of a channel with a fixed — i.e., indepen- 
dent of the input — pure output. We argued that the 
difficulty — that is, the low probability of success — in the 
realization of such a channel is essentially due to the ne- 
cessity of using for the output the same single photon by 
which the input logical state is encoded in the d modes. 
This is exactly because no photons are available in the 
ancillary ports. Obviously, if such a fixed pure output is 
readily available as ancillary state, the realization of the 
pure-fixed-output channel becomes trivial: the ancillary 
input state becomes the output. It is therefore evident 
that introducing non- vacuum ancillas would strongly af- 
fect the performance of our scheme. Another addition 
that we plan to consider is feed-forward, that is becom- 
ing a powerful and reliable tool in linear-optics quantum 
information processing [3, [3 ■ 

Another possible line of research is that of focusing on 
channels that are linked to interesting effects in quan- 
tum information processing. Indeed, our results can be 
thought of as a toolbox to be used in any optical exper- 
iment where some specific channel has to be applied, be 
it for the sake of simulating noise or for implementing a 
specific protocol. 

We expect our findings to trigger further theoretical 
studies on channel realization. In particular, we linked 
channel realization with more abstract notions of en- 
tanglement theory, and we hope that the study of less 
explored entanglement properties of states will conse- 
quently be stimulated. From a more practical point of 
view, our results provide a simple method for realizing 
arbitrary quantum channels using linear optics and stan- 
dard experimental techniques. Our results are ideal for 
experimental implementation relying on linear optics in 
combination with post-selection. While quantum chan- 
nels have been a widely discussed topic in theoretical 
quantum information, we expect our work to trigger an 
increased interest in the experimental study of this in- 
triguing topic. 

We thank Antonio Acin, Dominic Berry, Robert 
Prevedel, Kevin Resch and Armin Uhlmann for help- 
ful discussions. We acknowledge support from NSERC, 
Quantum Works, Ontario Centres of Excellence, CFI, On- 
tario MRI and ERA. 
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Appendix A: Proof of Lemma [T] 



Proof. By 
tion dni 



(PH]) in the defini- 
we obtam EciP) ~ 
1 - ^1 - C{i}},Y^ . Then, in or- 
der to obtain the lower bound it is sufficient to 



substituting 
_ of E<^{p) 



monotonically decreasing: 

E^{p) - max (l - - C{4>^f 



> max — I 1 — A / 1 — 




max V'piC(V'i 



(Al) 



The upper bound can be derived from the relation 
VI- x^>l-x. □ 

Appendix B: Proof of Theorem [3] 



Proof. We will use properties of the Shannon entropy 
Ei{{ri}) = — ^if'i\og2ri, defined for a probability dis- 
tribution {ri\. The von Neumann entropy of a quantum 
state p is equal to the Shannon entropy of its eigenvalues 
r^. In particular, for a pure bipartite state with Schmidt 
decomposition = X^i V^N)a|«)b, the entropy of 

the reduced one-party states p"^ and p^ is H{{\i}). For 
any pure ensemble {Pa,V'^^} we will denote by {A"}i 
the set of the squares of the Schmidt coefficients of tp-^^ , 
and define Aa,max = maxi{Af}i. For the entanglement of 
assistance it then holds 



= max VpaS'(pf) 
= max Vp,i/({A^},) 

, \ ^ TT ( \ ^ ^ Aq max 

S max > Patl A 

a. ma 

{Pa,^a}^ V 



1- A„ 



d- 1 



d-1 



< max hd 

{Pa,1>a} 



I / ^ J^a'^a, max I 
Id min y^PaXa.max 



a{A) 



(Bl) 



observe that vl 



is a concave function in x that is 



The first inequality is due to the fact that substitut- 
ing any subset of probabilities of some distribution with 
equally weighted probabilities can only increase the to- 
tal Shannon entropy. This is easily checked by know- 
ing that the flat probability distribution is the one with 
highest Shannon entropy, and that for any grouping of 
probabilities {r^} into two subsets {r,'^-*} and {r^''} of 
weight q and 1 ~ q, respectively, we have H{{ri}) = 
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hM + lH{{r^^^/q}) + (1 - q)H{{rf/{l - q)}). The 
second inequahty is due to the concavity of entropy. The 
second to last equality is due to the the monotonicity of 
hd in the interval [1/d, 1]. Indeed, X^aPa-^a.max > 
because Ao,max > l/d for all a. Finally, the last equal- 
ity comes from the relation (jllb| . Thus, using the fact 
that hd is invertible and monotonically decreasing in the 
range of interest, we obtain ct(A) < dh'^^(Eg,{J\)), i.e., 

ftucc > 1/ [dh^\E,{JA))]. 

For the upper bound we have 

EM = max Vp„iJ({An») 

{Pa,1pa} 

> max (-Pal0g2 (Aa,max)) 

{Pa,V'a} ^ 

- max (-logafVpaAa,, 
= -log2( min V'paAa.m 



(B2) 
(B3) 



los 



a(A) 



The first inequality comes from the fact that the min- 
entropy H^in{{ri}) = — log2 r-max of a probability distri- 
bution {rj, with r,„ax = max{ri}, satisfies H,^in{{ri}) < 
H{{ri}). The second inequality is due to the concavity 
of the logarithm. The second-to-last equality is due to 
the monotonicity of the logarithm. We finally arrive at 
the desired relation by exponentiation. □ 



Appendix C: Basis independence for the 
probabilistic constant-output channel 



Suppose t' = UtW; then 
J(A) 

= i^-p)Hd){^d\+P^' 

= {i-pM){i'+\+pUTU^ 



{U(g)U*) 



{i-pM){^i 



I 

2 

/ 

-pr (g) - 



{U(g)U 



where we have used the invariance of the maximally en- 
tangled state li^d) ~ , valid for all unitaries 
U. 



Appendix D: Decomposition saturating the bound 



One can write the Choi-Jamiolkowski state as the con- 
vex combination 



J(A) = (1 -p) 1^2+) (^2+ 1 + 2p(l ~s)-®- 

j ^ ^ (Dl) 
+ p(s~(l-s))|0)(0|®-, 
such that for the concave-roof of the geometric measure 
we find 



E2{J{A)>{l-p)E2m){^P+\) 



2p{l - s)E2 



+ p{s-{l~s))E2[\0){0\ 

= (l-p)i + 2p(l-s)i 
^ 1/2 -p{s- 1/2). 



Here we used that fact that £'g (|0) (0| ® //2) = 0— 
see the discussion just after Eq. (|18p — and that 
Eq{I/2®I/2) ~ 1/2, because the maximally mixed 
state of two qubits can be seen as the convex combination 
of pure maximally entangled states. 



